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SUMMARY 


A  theoreticel  and  experimental  investigation  was  undertaken 
of  heat  transfer  effects  at  the  surface  of  a  flat  plate  on  which 
combustion  is  occurring  in  a  supersonic  air  stream.  A  simplified 
theoretical  model  is  developed  which  leads  to  expressions  for  the 
heat  transfer  coefficient  and  heat  transfer  rate  at  the  surface. 

The  experimental  technique  used  to  obtain  heat  transfer  parameters 
is  described.  The  experimental  tests  were  conducted  at  Mach  number 
3.1  using  the  pyrophoric  fuel  triethylalu ainum.  A  co.aparison 
between  the  theoretical  and  experimental  results  indicates  that  the 
former  yields  or;  'r-of-magnitude  estimates  of  the  heat  transfer 
coefficient  and  rate  whicli  are  satisfactory  for  purposes  of 
engineering  design.  The  theory  is  then  used  to  obtain  estimates  of 
heat  transfer  parameters  at  the  surface  of  a  cone-shaped  interceptor 
executing  a  coast  trajectory.  The  results  are  shown  to  be  well 
correlated  by  an  empirical  expression  similar  to  that  encountered 
in  the  literature  for  turbulent  non-reacting  flow  over  flat  plates 
having  both  constant  and  variable  surface  temperatures. 
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I.  INTRODUCTION 

A  number  of  schemes  exist  for  utilizing  a  combustion  process 

1  2  3-8 

to  produce  aerodynamic  lit l.  forces  *  ,  vehicle  control  forces  ^ 

9-11 

and  thrust  vector  control  of  rocket  motors  .  In  each  case, 
the  underlying  principle  is  the  same:  To  inject  and  burn  a 
combustible  liquid  or  gas  at  some  surface  in  a  supersonic  flow 
environment,  leading  to  a  deflection  of  sti'eamlines  and  attendant 
pressure  increases. 

One  of  the  problems  introduced  by  the  application  of  combustion 
at  a  surface  is  the  need  for  a  priori  estimates  of  the  surface 
heat  transfer  to  ensure  an  appropriate  design.  In  the  following 
section,  an  approximate  theoretical  analysis  is  developed  for 
calculating  the  heat  transfer  to  a  turbulent  flat  plate  in  the 
presence  of  heat  generated  by  combustion.  In  an  "-ffort  to 
detennine  the  effectiveness  of  this  analysis  in  predicting  heat 
transfer  rates  and  coefficients,  a  series  of  experiments  was 
undertaken  in  which  the  pyrophoric  fuel  triethylaluminum  (TEA) 
was  injected  and  burned  at  the  surface  of  a  flat  plate  in  Mach  3.1 
flow.  The  details  of  these  experiments  are  contained  in  Reference 
4,  and  are  briefly  s\immarized  in  Section  III. A  of  this  report. 

During  these  tests,  the  plate  surface  temperature  was 
monitored  by  a  heat  transfer  gage  of  special  design,  which  is 
described  in  Section  III. A.  From  these  readings,  heat  transfer 
coefficients  and  heat  transfer  rates  at  the  plate  surface  were 
deduced  (Section  III.B) .  The  experimental  results  obtained  by 
this  means  are  presented  in  Section  III.C. 

The  theoretical  and  experimental  results  are  compared  in 
Section  IV.  As  an  example  of  the  use  of  the  theoretical  analysis 
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in  obtaining  estimates  of  heat  transfer  coefficients  and  heat 

transfer  rates,  these  parameters  are  calculated  for  the  case  of 

external  burning  employed  to  generate  a  control  force  on  a 

cone-shaped  interceptor  executing  a  coast  trajectory  (Section  V) . 

Throughout  owcse  calculations,  the  fuel  considered  is  TEA? 

however,  this  fuel  is  sufficiently  similar  to  most  hydrocarbon 

fuels  (H  fid  18,000  BTU/lbm)  to  permit  the  use  of  these  results 
c 

in  estimating  heat  transfer  coefficients  and  rates  on  vehicles 
executing  similar  maneuvers  and  using  other  hydrocarbon-like 
fuels . 


II.  THEORETICAL  HFAT  TRANSFER  ANALYSIS 


TR-486 

3 


An  approximate  model  of  heat  transfer  to  a  flat  plate  in 

the  presence  of  combustion  may  be  formulated  by  a  modification 

12 

of  the  analysis  developed  by  Mayer  .  The  two-dimensional  model 
employed  is  depicted  in  Figure  1.  Here,  conditions  denoted  by 
the  subscript  co  refer  to  conditions  in  the  flame  zone  where 
the  stagnation  temperature  is  varying  with  x,  rather  than  the 
inviscid  nozzle  flow  field  examined  by  Mayer.  Effects  such  as 
the  deposition  of  insulating  solids  on  the  surface  and  radiation 
are  neglected. 

Referring  to  Figure  lb,  the  various  energy  terms  ent  '.ring 
the  control  volume  are  understood  to  include  both  sensible  and 
chemical  energies.  Thus, 


c  T  dy  +  dh 
p  o 


c,l 


(1) 


where  h  is  the  total  sensible  enthalpy  (including  kinetic 
s 

energy)  and  h^  is  the  chemical  enthalpy.  In  like  manner, 


=c  T 
p  c,  or; 


_d 

dx 


P  I 

j 


dx  +  dh  - 

c, 


(2) 


uc  T  dy  + 
P  o 


O  uc  T  dy 
po 


dx 


(3) 


+  dh  _ 
c,  3 

Application  of  the  conservacion  of  energy  to  the  control 
volume  leads  to 


q  =c  T 
^  p  o, 00 


r  T  1 

d  r 

,  o 

dx  7  ^  ^ 

"  T 

~'o 

O,  CB 

dy  + 


dx 


(4) 
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where  dq^  represents  the  net  amount  of  energy  entering  the 
control  volume  by  virtue  of  the  addition  of  chemical  energy. 
At  the  edge  of  the  boundary  layer,  the  gas  dynamic  conditions 
approach  conditions  in  the  flame  zone.  If  it  is  assumed  that 
the  flame  zone  extends  well  above  the  plate  surface,  and  that 
conditions  in  this  region  are  independent  of  y.  Equation  (4) 
may  oe  written  to  a  fair  degree  of  approximation: 


A  "stagnation  temperature  thickness"  is  introduced 
defined  by 


which  has  dimensions  of  length.  The  thermal  Reynolds  number 
is  then  defined  as 


P  u 

'^OO  OD 


r 


(7) 


Introducing  these  terms  into  Equation  (5),  and  initially 
assiaming  that  c^,  \i,  and  Pr  are  constants  at  flame  zone  conditions 
(the  case  of  variable  /i,  and  c^  and  Pr  taking  on  constant 
values  other  than  oo  will  be  treated  subsequently) ,  leads  to 

d 

^  oo  p  0,00  dx  r  dx  '  ^ 

oo 
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Introducing  the  adiabatic  heat  transfer  coefficient  and  Stanton 
number,  respectively: 


h 


T 


ad  w 


(9) 


St  = 


0  c 

'^QO  p  U 

00  00 


(10) 


Equation  (8)  becomes 


P  U 

d  „  *^00  00 

Re_  =  - 

dx  r  H 

'^ao 


T  ,-T 
ad  w 

T 

0,00 


St  -  q* 
c 


(11) 


in  which 


dq 


q*  =  - 

c  ^  c 


OD  p  T 

00  0,00 


c 

dx 


(12) 


as 


A  semi-empirical  recovery  factor  is  now  employed,  defined 


R  s 


T  ,  -  T 
ad  oo 

T  -T 
0,00  OD 


(13) 


where  R  -  for  uurbulent  boundary  layers^^.  Using  the 

relation 


T  =  T 

0,00  oo 


1  + 


2zi 

2 


(14) 


in  Equation  (13)  leads  to  the  foil  ;wing  expression  for  the 
adiabatic  wall  temperature: 


rp  — m 

"ad  "^0,00 


V— 1  2 

^  ^  0  Mj+i 

_ _ <p 

V'-l  2 
1+  M 
2  oo 


(15) 


Following  Mayer,  we  now  assume  that  the  relation  between  Re 
and  St,  which  is  necessary  to  solve  Equation  (11),  is  the  same 
relation  which  may  be  derived  for  the  case  of  uniform  flow  along 


a  flat  plate.  This  latter  expression  may  be  obtained  as 
follows : 

The  Blasius  equation  for  the  skin  friction  on  a  flat  plate 
for  constant  fluid  properties  is 


f#  oo  _  _  -b 
—A —  =  B  Re 
2  X,GD 


(16) 


A  modified  Reynolds  analogy  is  employed: 


f.oo  -2/3  -2/3  -b 

St  =  -■  i'  Pr  ^  =  B  Pr  Re 

OO  2  OO  00  X,(D 


(17) 


where  the  semi-empirical  Blasius  parameters  have  the  values 
b  =  1/5,  B  =  0.0296  for  a  turbulent  flat  plate  boundary  layer, 
and  the  Reynolds  number  is  given  by  Re  =  o  U  x/u 

X ,  CD  OD  OO  CD 

Substitution  of  Equation  (17)  into  Equation  (11)  gives 


r, 

dx  r 


P  u 

'^CD  CD 


T  ,-T 
ad  w 


0,CD 


BPr”^^^  Re"^  -q*  (18) 

X,(X)  c 


In  order  to  integrate  Equatidn  (18),  the  following 
simplifying  assumptions  will  be  made: 

1.  The  term  (T  j-T  )/T  will  be  assumed  to  be  independent 

aa  w  0,00  ^ 

of  X.  This  assumption  was  employed  by  Mayer  as  well  as  for  flow 
through  a  rocket  nozzle,  and  an  a  posteriori  examination  of  the 
results  generated  by  this  analysis  indicates  that  this  parameter 
changes  by  only  4  or  5%  in  going  from  (x/L)=0  to  (x/L)=l,  where 
L  is  the  plate  length. 

2.  In  integrating  the  second  term  in  Equation  (18),  T 

o,CD 

will  be  assumed  independent  of  x.  The  justification  for  this 
assumption  is  that  a  comparison  between  the  analytical  and 
experimental  results  (Section  IV)  indicates  that  the  local 
equivalence  ratio  (fuel-air  ratio/fuel-air  ratio  at  stoichiometric) 


I 

I 

I 

1 
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in  the  vicinity  of  the  plate  is  substantially  less  then 


unity,  resulting  in  relatively  small  increases  in  T  as  x 

0,00 

increases.  This  also  has  the  advantage  of  delaying  the  necessity 
for  specifying  q^(x),  thereby  generalizing  the  final  result. 

Employing  these  two  assumptions,  integration  of  Equation  (18) 
yields 


0,00/ 


T,  -2/3  „  (1-b) 
B  Pr  Re 

CD _ X,  oo 

1-b 


U  c  T 

P  O  ,  00 
^OO 


(19) 


where  the  initial  conditions  that  Re_  and  q  are  zero  at  x=0 

r  c 

were  used. 

Using  Equation  (17)  in  Ecjuation  (19)  leads  to 


ROj,  +Ch  = 


_ 00 _ 

1-b 


IT  -T  ' 
ad  w 


CO 


(1-b) 

(St  )  -b 
OD 


(20) 


in  which  the  dimensionless  "chemical  energy  number"  is  given  by 


Ch  s 


a  c  T 

CD  p  o ,  00 

^oo 


(21) 


Equation  (20)  may  now  be  used  in  Equation  (11)  to  obtain 
the  following  first-order  differential  ecpiation: 


(1-b) 

-b 


P  u 

'^CD  00 


CD 


iizM. 


/T,  r,  -2/3, 1/b 
(B  Pr  ) 

CO 


(22) 


which,  upon  integration  and  using  Equation  (10) ,  leads  to 


h  = 


B  Pr"^^^  P  U 

CD  '^OO  (D  CD 


(23) 


Re 


X,  00 
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The  initial  condition  St - s*<B)(Rej=0)  at  x=0  was  employed  above. 

The  effect  of  variable  density  and  viscosity  may  now  be 
introduced  in  the  following  manner: 

We  assume  that  all  fluids  are  ideal  gases,  so  that 

fi  «  t"^  (24) 

and  that  the  viscosity  is  a  function  of  temperature  only;  this 
functional  relationship  may  be  written 

tl  ^  (25) 

where  to  is  an  empirically-determined  parameter.  The  Prandtl 

number  and  specific  heat  are  assumed  to  be  constant  at  some 

mean  value  between  T  and  T  (i.e.,  Pr=Pr,c  =c  ). 

0,00  w  P  P 

We  now  employ  Eckert ' s  reference  temperature  method  as 

12 

suggested  by  Mayer  ,  in  which  the  modified  Stanton  number  (to 
account  for  the  variable  fluid  properties)  may  be  written 


1-b 

) 


(26) 


where  /i*  and  P*  are  evaluated  at  an  empirical  reference 
temperature  given  by 


T*  =  ^2  (T  +T  )  +  0.22  (T  ^-T  ) 
w  CD  ad  (a> 


(27) 


Using  Equations  (24)  and  (25)  in  Equation  (26)  yields 
St 


St 


T  l-b(l+ai) 

=  (  —-  )  =  i8 


CD 


T* 


(28) 
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Therefore,  Equation  (17)  and  (28)  lead  ^o 


St  =  St 


^  BVt 


-2/3 


CD 


CD 


(29) 


Re 


'X,  CD 


Substituting  Ecjuation  (29)  into  Equation  (11)  and  integrating 
the  resulting  expression  gives 


„  ,  ad  w 

Rep  +  Ch  =  ^ - 

0,00  I 


O  B  7r"^^^) 
CD 


1-b 


lA 


St 


-b 


(30) 


Where  Equation  (17)  was  also  used.  Once  again,  the  initial 

condition  Re^,  =  Ch=0  at  x=0  was  employed. 

Resubstituting  this  result  back  in  Equation  (11)  gives 

the  following  first-order  differential  equation: 

,  ^  (1-b) 

T  ,  lA  -7“  ,  ,  p  U 

J=.  .  Q  St  =  - — -  •  pnN 

St  dx  ^  --2/3,1A  II 


(B  Pr  '  ) 


CD 


which  may  be  integrated  (employing  a  substitution  of  variables. 


G=/3  St)  to  give 


h  = 


-2/3 

a  1-b  ~  ,  p  C  U  . 

^  J _ (B  Pr  )  (  ^co  p  CD  ) 


/ 


1  tt  -1  ^ 

P  U  U  dx 
1-b  ^CD  CD  '^CD 


(32) 


where  Equation  (10)  and  the  initial  condition  St- — ••cd  at  x=0 


were  used. 
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It  is  of  interest  to  note  that  Equations  (23)  and  (32)  are 
identical  in  form  to  the  equivalent  expressions  derived  by 
Mayer.  This  stems  from  the  two  m.ajor  assumptions  mentioned  above 
which  permitted  the  lumping  of  the  thermal  Reynolds  number, Re  , 
and  the  chemical  energy  number,  Ch,  into  an  effective  single 
dimensionless  grouping,  (R^P  +Ch) .  For  the  case  of  a  flat  plate, 
and  substituting  the  empirically-determined  values  of  b  and  B, 
Equation  (32)  becomes 


h 


0.0296  jS  p  c  U 

_ _ _ *^CP  p  OP 


Re 


1/5 

X,  00 


(33) 


and  the  heat  transfer  rate  to  the  wall  may  be  obtained  from 
Equation  (9) .  It  remains  to  specify  some  suitable  relation  for 
^  (x) ,  which  may,  of  course,  be  empirically-determined, 
and  to  choose  the  viscosity  relationship,  Eq\  ation  (25). 

For  simplicity,  we  assume  a  relationship  between  the 
chemical  heat  input  and  the  flame  zone  stagnation  temperature  of 
the  form 


xq 


c 


c 

P 


(T 


o,  oo 


which  may  be  manipulated  to  yield 


(34) 


T 

0,00 


rifli 

c 


c 

p 


+  T 

0,1 


(35) 


where  fj  is  the  equivalence  ratio,  f  is  the  stoichiometric  fuel- 

air  ratio,  and  H  is  the  heat  of  combustion  of  the  fuel. 

c 
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A  reasonable  empirical  expression  for  the  viscosity- 
temperature  relationship  is 

*7  r\  ^  Cs. 

li  =  0.835  X  10“  •  T  Ibm/lt-sec  (36) 

00 

which  is  shown  in  Figure  2  along  with  viscosities  of  the  species 
likely  i o  be  found  in  large  quantities  in  the  flame  zone.  It 
may  be  seen  that  Equation  (36)  provides  a  suitable  average 
flame  zone  viscosity  over  a  wide  range  of  temperature. 

It  should  be  noted  that  in  the  above  analysis  thermo- 

13 

chemical  and  kinetic  effects,  of  the  type  examined  by  Rosner 

have  not  been  included.  It  will  be  seen  (Section  IV)  that  the 

simplified  model  presented  above  yields  results  which  are  in 

reasonably  good  quantitative  agreement  with  experimental  results, 

and,  as  a  consequence,  supplies  estimates  of  the  heat  transfer 

coefficient  and  heat  transfer  rate  which  are  sufficiently  accurate 

for  most  engineering  purposes. 

It  will  be  shown  in  Section  IV  that  for  given  flame  zone 

conditions,  neither  the  equivalence  ratio,  molecular  weight, 

nor  Prandtl  number  substantially  affect  the  value  of  heat 

transfer  coefficient,  and  that  only  7]  affects  the  heat  transfer 

rate  to  the  wall  to  a  substantial  extent.  This  latter  effect  stems 

from  the  direct  relationship  between  T  and  77  through  Equations 

ad 

(15)  and  (35).  The  insensitivity  of  h  to  the  parameters  mentioned 
above  permits  a  wide  latitude  in  the  estimation  of  these 
parameters  and  circumvents  the  need  for  an  extensive  chemical 
equilibrium  analysis  of  the  combustion  process  to  determine 
accurate  values  for  9?Z  a^^-d  Pr  . 
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III.  EXPERIMENTAL  TECHNIQUES 

A.  Description  of  Experimental  Apparatus 

The  experimental  program  was  conducted  in  the  GASL 
Mach  number  3.1  high-temperature  blowdown  wind  tunnel  using  the 
pyrophoric  fuel  triethylaluminum  (TEA) .  The  nominal  operating 
conditions  were  100-250  psia  stagnation  pressure  and  1500-2500°R 
stagnation  temperature.  The  tunnel  is  described  in  greater  detail 
in  Reference  4.  TEA  was  injected  through  an  axially-aligned  , 
triangular- shaped  pattern  of  injector  holes  shown  with  the 
model  in  Figure  3. 

Since  the  stagnation  conditions  were  not  invariant 

with  time,  a  transient  heat  transfer  measuring  technique  was 

employed.  The  conventional  thin-walled  model  technique  which  is 

frequently  used  in  connection  with  transient  heat  transfer 

experiments  is  not  applicable  here  because  of  wind  tunnel 

limitations.  The  time  for  the  tunnel  to  achieve  steady  flow  is 

approximately  five  seconds;  therefore,  it  is  impossible  to 

obtain  a  temperature  history  which  reflects  only  the  steady- 

flow  conditions.  Thus,  for  the  heating  rates  available  the 

model  must  be  thick  and  net  thin.  It  was  therefore  decided 

14 

to  use  the  method  discussed  by  Crep::i  and  Libby  .  This 
technique  is  also  a  transient  one  and  involves  a  cylindrical 
plug  attached  to  the  model  at  the  surface  but  insulated  from 
it  for  the  rem.aining  portion  of  its  length  by  an  air  gap  (Figure  4)  . 
It  is  assumed  that  this  gap  acts  as  a  sufficiently  good  in¬ 
sulator  to  prevent  any  heat  transfer  between  the  plug  and  the 
surrounding  body,  thus  permitting  heat  flow  only  parallel  to 
the  plag  axis.  This  arrangement  transforms  the  unsteady  heat 
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conduction  problem  from  one  involving  axial  and  radial 
conduction  into  one  which  is  locally  only  one-dimensional 
and  analogous  to  a  rod  heated  at  one  end.  In  general,  if  the 
surface  temperature  history  and  the  initial  temperature  distri¬ 
bution  in  the  plug  are  known,  the  solution  for  the  rate  of  heat 
transfer  to  the  plug  can  be  determined.  The  initial  temperature 
distribution  in  the  plug  is  assumed  to  be  constant  and  care  is 
taken  to  establish  such  a  state  prior  to  testing.  The  surface 
temperature  history  is  determined  by  a  thermocouple  located 
at  the  plug  surface. 

The  junccion  of  the  thermocouple  is  made  as  close 
to  the  surface  as  possible  by  the  following  technique:  The 
individual  thermocouple  wire  insulation  is  drawn  up  through 
a  small  hole  in  the  plug  at  the  surface,  then  a  tapered  pin 
of  the  same  material  as  the  model  and  plug  is  driven  into  the 
hole  sc  as  to  cut  the  insulation  against  the  wall  of  the  plug 
at  the  surface.  The  pin  is  filed  off  and  the  surface  smoothed. 
The  therTTivOcouples  are  then  checked  to  detect  buried  junctions 
by  applying  a  small  heat  input  and  observing  the  response  time. 

B.  Heat  Transfer  Gage  Analysis 

For  relatively  low  rates  of  heat  transfer  or  short 
periods  of  heating,  the  thermocouple  plug  can  be  considered  a 
semi-infinite  solid.  That  is,  the  inside  surface  of  the  plug  is 
assumed  to  remain  at  essentially  the  same  temperature  as  it  was 
at  the  initiation  of  heating.  To  obtain  a  general  solution  to 
the  heat  conduction  equation ^ temperature  at  the  surface  is 
considered  an  arbitrary  function  of  time.  This  method  gives 
local  values  at  each  interval  of  time  and  eliminates  the  effect 
of  non-uniform  starting  conditions. 
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For  the  following  boundary  conditions : 

T  =  T.  at  t=0 
1 

T  “  T.  +  (p  (t)  at  y=0 

T  =  T.  at  y=0D 
1 

14 

it  can  be  shown  that  by  integrating  the  partial  time  derivative 
of  the  temperature  from  zero  to  infinity,  the  heat  transfer  is; 


S<=o 


(t-X)"^^^  dX 


(37) 


For  numerical  integration  and  data  reduction.  Equation  (37) 
^an  be  approximated  by 


2^J1"  (p'  (t) 


^  <p(t-0 

y/2 


where  (>0  and  is  selected  such  that 


t-c 

f 

j 

o 


2(t-X) 


3/2 


dX 


(38) 


(p'  (X)  -  <p'  (t)  for  t-c  ^  X  ^  t 

Thus,  for  any  measured  surface  temperature  history,  the  aero¬ 
dynamic  heat,  transfer  rate  can  be  determined  from  Equation  (38) . 
C.  Experimental  Results 

Tests  were  conducted  at  Mach  number  3.1  and  nominal 
stagnation  pressures  of  IIC  and  250  psia.  The  basic  data  were: 
heat  transfer  coefficient,  heat  transfer  rate,  surface  tempera¬ 
ture  distribution,  stagnation  temperature,  and  injectant  mass 
flow  rate.  The  variation  of  these  parameters  with  time  is  pre¬ 
sented  'u  Figures  5  and  6. 
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It  is  observed  in  Figure  5  that  for  tests  conducted 

at  ilO  psia  stagnation  pressure  (p.=2.63  psia)  and  stagnation 

o  ^ 

temperature  approximately  2300  R,  the  centerline  heat  transfer 
coefficient  and  heat  transfer  rate  2.5  inches  downstream  of  the 
injector  stabilij'.e  at  some  value  and  then  decrease  with  increased 
injectant  mass  flow.  This  point  is  significant  since  it  appears 
that  for  given  stagnation  conditions  there  exists  an  optimum 
injectant  mass  flow  rate  which  if  exceeded  results  in  a  decrease 
in  h  (x)  and  q^(x)  .  A  similar  existence  of  ■^n  optimum  mass  flow 
was  noted  earlier  in  the  force  productic  i  mechanism,  described 
in  Reference  3,  and  was  attributed  to  cooling  of  the  conioustion 
gases  by  unburned  fuel.  This  cooling  coupled  with  the  corres¬ 
ponding  local  pressure  decrease  due  to  droplet  breakup,  e'/apora- 
tion  effects  and  local  aerodynamic  expansions  tend  to  lower 
the  efficiency  of  the  combustion  process  and  result  in  a 
corresp'^nding  decrease  in  h(x)  and  q^(x),  shown  as  functions 
of  injectant  mass  flow  in  Figure  7a. 

Tests  at  200  psia  stagnation  pressure  (p^=4.7  psia) 
and  1500-1900°R  stagnation  temperature  resulted  in  opposite 
trends  in  the  heat  transfer  ra  i  and  heat  transfer  coefficient 
data  6.5  inches  downstream  of  the  injector  (Figure  6  and  7b). 

This  result,  which  was  partially  due  to  decreasing  stagnation 
temperature,  demonstrates  the  greater  dependence  of  g^(-) 
ignition  delay  and  chemical  reaction  times.  This  effect  is 
explained  in  detail  in  Section  IV. 

The  development  of  the  ccir:bustion  process  along  the 
surface  of  the  model  may  be  traced  by  reference  to  the  axial 
gradient  of  h (x)  and  (Figure  8).  For  the  low  pressure 

case  this  gradient  initially  decreases,  indicating  a  constant 
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or  decreasing  effective  local  equivalence  ra^^io,  and  then  in¬ 
creases  do^vnstream  as  the  fuel  spreads  on  the  surface  of  the 
model  and  mere  air  becomes  available.  In  the  high  pressure 
case  the  gradient  continuously  increases,  indicating  a 
steadily- increasing  local  equivalence  ratio. 

Finally,  in  Figure  9  it  is  shown  that  the  reaction 
is  approximately  symmetric  about  the  centerline  and  also  that 
a  more  efficie.  t  burning  process  occurs  as  one  proceeds  from 
the  centerline  toward  the  outer  edge  of  the  flat  plate  model, 
(in  the  directic/i  of  increasing  z) .  This  latter  phenomenon 
implies  that  a  fuel-rich  condition  most  likely  exists  along  the 
plate  centerline,  which  becomes  more  nearly  sto ! chiomerric  in 
proceeding  to  the  outer  boundary  of  the  flame  zone. 
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IV.  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT 

In  order  to  determine  the  effectiveness  of  the  theoretical 

analysis  (detailed  in  Section  II)  in  predicting  heat  transfer 

coefficients  and  rates,  a  comparison  for  two  typical  experimental 

runs  is  made  in  this  section.  The  test  conditions  chosen  for 

comparison  are  tabulated  in  Figure  10-13.  The  times  were  chosen 

as  being  representative  of  periods  when  a  quasi-steady-state 

was  achieved.  (Section  III,  Figures  5  and  6). 

Comparison  between  the  theoretical  and  experimental  Vclues 

of  h(x)  and  two  typical  runs  are  shown  in  Figures 

10  through  13.  It  may  be  seen  that,  in  spite  of  the  limitations 

imposed  by  the  assumptions  used  in  the  analysis  of  Section  II, 

both  h (x)  anc  q  (x)  agree  with  the  experimental  values  to  with- 

w 

in  an  order-of-magnitude  when  a  fixed  value  of  rj  =  0.01  is  used, 

in  spite  of  the  fact  that  the  experimental  flame  was  not  two- 
4 

dimensional  .  It  will  be  seen  presently  that  the  theoretical 

value  of  h(x)  is  fairly  insensitive  to  rj.  However,  q  (x)  depends 

upon  77  to  a  much  larger  extent.  When  a  low  value  of  t)  is  chosen 

77  =  O.Ol,  the  values  (Figures  11  and  13)  are  in  reasonably 

good  quantitative  agreement.  The  corresponding  values  of  h(x) 

for  this  value  of  77  are  shown  in  Figures  10  and  12. 

The  differences  in  the  sign  of  the  slopes  of  the  theoretical 

and  experimental  curves  most  likely  stems  from  a  variation  in 

the  value  of  the  equivalence  ratio  with  plate  length  in  the 

3  4 

experimental  case.  That  is,  there  is  experimental  evidence  ' 
which  indicates  that  the  overall  chemical  reaction  accelerates 
as  X  increases,  because  of  droplet  breakup,  evaporation  effects, 
and  changes  in  the  ignition  delay  and  chemical  reaction  times 
resulting  from  variations  in  the  local  environmental  conditions. 
The  variable  77  curve  in  Figure  13,  in  which  a  linear  increase  of 
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71  with  X  was  arbitrarily  taken,  demonstrates  that  theoretical 
curves  with  slopes  having  the  same  sign  as  the  experimental 
curves  can  be  generated  by  including  in  the  calculation  this 
additional  experimentally-observed  effect.  However,  the  order 
of-magnitude  agreement  obtained  with  constant  77  is  regarded 
as  sufficiently  accurate  to  provide  satisfactory  engineering 
estimates  of  heat  transfer  coefficients  and  rates. 

In  obtaining  the  theoretical  estimates  of  h  (x)  and 
shown  in  Figures  10  through  13,  the  follov/ing  values  of  input 
parameters,  characteristic  of  TEA  combustion  and  the  experi¬ 
mental  conditions,  were  emp‘'oyed: 

fmr.*  =  0.08  Ibm./lbm 
TEA  f  a 

H  =  18,300  BTU/lbm^ 

c,TEA  f 

Pr  =  0.75 

c  =0.30  BTU/lbm-°R 
P 

^  =  28  Ibm/lb-mole 

Tj  =  0.01 

T  .  =  2300°R 
0,1 

p  =  4.2  psia 

<30 

M  =  1.25 
00 

L  =  10.25  in. 
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The  chosen  value  of  represents  the  average  of  the  values 

at  the  plate  surface  as  measured  by  a  large  nxamber  of  static 

pressure  taps.  The  value  of  was  chosen  as  follows: 

In  these  experiments,  the  combustion  process  forms  an 

approximate  half-cone  region  which  tends  to  divert  the  stream- 

4 

lines  in  much  the  same  way  as  a  solid  half-cone  .  By  measuring 

14 

the  cone  angle  from  photographs  and  using  cone  charts  ,  the 

approj^imate  Mach  number  at  the  edge  of  the  burning  zone  can  be 

5 

calculated.  Billig  shows  that  the  variation  of  Mach  number 
from  the  wall  through  the  combustion  region  is  approximately 
linear.  As  a  result,  an  average  Mach  number  in  this  zone  may 
be  taken  as  one-half  the  value  ac  the  outer  edge  of  the  burning 
region. 

Clearly,  further  manipulation  of  such  parameters  as 

and  M  can  bring  the  theoretical  and  experimental  values  of  h(x) 

CD 

and  still  closer.  However,  for  present  purposes,  such 

empirical  manipulation  is  not  deemed  warranted. 

It  is  instructive  to  examine  the  effect  of  certain  para¬ 
meters  on  the  values  of  h  and  q^.  Figure  14  shows  the  effect  of 
varying  the  equivalezice  ratio,  rj.  As  was  previously  mentioned, 
changing  t)  by  a  factor  of  ten  results  in  a  relatively  small 
change  in  h;  however,  the  wall  heat  transfer  varies  substantially. 
As  pointed  out  in  Section  II,  this  latter  effect  stems  from  the 
dependence  of  on  77  through  Equations  (15)  and  (35) . 

The  effect  of  a  specific  heat  variation  is  shown  in  Figure  15. 
As  may  h<r  seen,  a  20%  variation  in  c  leads  to  a  corresponding 
change  in  h  and  q^,  primarily  through  the  direct  dependence  on  h 

on  c  in  Equation  (33).  However,  the  flame  zone  specific  heat 
P 

for  hydrocarbon- like  fuels  changes  relatively  little  over  wide 
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ranges  of  pressure,  temperature  and  equivalence  ratio  .  As 
a  result,  this  parameter  can  be  estimated  with  reasonable 
accuracy  from  existing  data. 

The  effects  of  molecular  weight  and  Prandtl  number  changes 
are  shown  in  Figures  16  and  17.  In  both  cases,  h  and  change 
very  little,  indicating  that  highly  accurate  estimates  of  these 
parameters  are  not  required  for  present  purposes. 

The  fai.rly  wide  margin  of  uncertainty  which  can  be  toleifated 
in  the  estimation  of  input  parameters,  adds  to  the  utility  of  the 
theoretical  analysis  in  estimating  heat  transfer  coefficients 
and  rates  encountered  on  the  surface  of  a  flight  vehicle  on  which 
combustion  is  occurring.  Such  estimates  are  discussed  in  the 
next  section. 
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V.  THEORETICAL  HEAT  TRANSFER  RESULTS  CORRESPONDING  TO 

FLIGHT  CONDITIONS 

The  mission  selected  for  evaluation  involves  the  uniformly- 
decelerating  coast  trajectory  executed  by  an  8°  half-angle  cone 
from  10,000  fps  at  3000  foot  altitude  to  6000  fps  at  20,000  feet. 
The  burning  at  the  cone  surface  is  assumed  to  produce  an 
approximate  half-cone  of  15°  included  angle,  from  which  the 
average  flame  zone  Mach  number  may  be  calculated  as  in  the 
previous  section.  The  surface  pressure  in  the  disturbed  region 
has  been  taken  as  twice  the  undisturbed  cone  surface  pressure, 
and  the  remaining  input  parameters  were  chosen  as  in  the  previous 
section.  The  characteristic  length  of  the  system,  L,  was 
taken  as  five  feet. 

The  results  are  shown  in  Figures  18  and  19.  Once  again, 
the  lack  of  sensitivity  of  the  heat  transfer  coefficient  to 
the  input  parameters  {r}  and  T^  in  this  case)  is  demonstrated  in 
Figure  18,  where  h  is  seen  to  fall  within  a  relatively  narrow 
band  for  the  ranges 

500°R  ST  S  4000°R 
w 

0.01  s  7]  s  1.00 

For  computational  purposes  involving  approximately  similar 
trajectories  the  heat  transfer  coefficient  may  be  selected  from 
Figure  18. 

The  corresponding  values  of  heat  transfer  rates  to  the 
vehicle  wall  are  shown  in  Figure  19.  The  greater  sensitivity 
of  to  the  input  data  may  again  be  seen.  However,  the  presented 
data  may  be  used  as  an  order-of-magnitude  estimate  of  the  wall 
heat  transfer  rates. 
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It  should  be  noted  in  Figure  19  that  for  high  values  of 
\  and  low  values  of  r?  at  high  altitude  the  heat  transfer  to 
the  wall  becomes  negative.  This  corresponds  to  the  case  where 
"^ad  ^  "^w  the  case  where  the  wall  has  an  internal  heat- 

generation  mechanism.  Clearly,  in  the  present  context  this  is 
not  a  physical  solution. 

In  the  utilization  of  heat  transfer  data  the  dimensionless 
Nusselt  number  is  often  used  in  the  presentation  of  results. 
This  parameter  is  defined  as 


where  k  is  the  thermal  conductivity  of  the  gas.  In  evaluating 
the  thermal  conductivity,  the  Eucken  expression  for  polv- 
atomic  gases  at  low  density  was  used^^: 


k  "  (c_  +  -r 


R 


4  55^ 


)  M 


where  the  viscosity,  ^  ,  was  computed  from  Equation  (36).  As 
may  be^seen  in  Figure  20,  the  Nusselt  number  varies  from  about 
1  X  10  to  4  X  10  for  the  trajectory  under  consideration. 
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VI .  CONCLUDING  REMARKS 

An  approximate  analysis  of  heat  transfer  to  a  flat  plate 
in  the  presence  of  combustion  has  been  formulated.  The  results 
generated  by  this  analysis  are  in  reasonable  quantitative  agree¬ 
ment  with  the  results  of  a  series  of  experiments  conducted  in  a 
Mach  3.1  nozzle, using  the  pyrophoric  fuel  triethylaluminum. 

Although  it  was  initially  assumed  that  a  constant  value  of 
equivalence  ratio  exists  along  the  plate  surface,  it  was  found 
that  the  trend  of  the  theoretical  curves  best  correspond  with 
the  experimental  trends  when  the  analysis  is  modified  to  include 
th''  effect  of  variable  Tj(drj/dx  >0).  In  any  event,  it  may  be 
concluded  on  the  basis  of  this  comparison  that  the  analysis 
will  yield  heat  transfer  ra^-es  and  heat  transfer  coefficients 
which  are  in  order-of-magnitude  agreement  with  the  experimental 
values,  and  are  sufficiently  accurate  for  most  engineering  design 
calculations. 

It  was  found  that  the  theoretical  values  of  h  and  q^  are 
fairly  insensitive  to  the  local  values  of  molecular  weight  and 
Prandtl  number.  In  addition,  the  heat  transfer  coefficient  was 
found  to  be  only  a  weak  function  of  the  equivalence  ratio,  wh  le 
q^  depends  upon  t)  to  a  much  larger  extent,  through  the  dependence 
of  the  adiabatic  wall  temperature  on  the  progress  of  the  combustion 
process.  Although  both  heat  transfer  parameters  are  strong 
functions  of  the  local  specific  heat,  this  latter  parameter  varies 
over  a  relatively  restricted  range  and  may  be  estimated  with 
sufficient  accuracy  by  reference  to  the  literature  (e.g.. 

Reference  15) . 
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A  relatively  simple  experimental  technique  for  obtaining 
values  of  heat  transfer  coef ficienctj  and  rates,  and  the 
theoretical  analysis  which  suggests  this  technique  were  discussed. 
The  experimental  results  serve  as  a  basis  for  the  above  evalua¬ 
tion  of  the  theoretical  analysis  developed  in  Section  II. 

Finally,  estimates  were  made  of  h  and  q^  at  the  surface  of 
a  cone-shaped  interceptor  executing  a  coast  trajectory  in  the 
altitude  range  3000 ' -20, COO ' .  It  was  assumed  that  an  external 
burning-control  force  generation  process  is  occurring  at  the 
cone  surface,  and  that  the  vehicle  is  sufficiently  large  that 
local  curvature  effects  may  be  ignored.  It  was  found  that  for 
extensive  ranges  of  equivalence  ratio  and  wall  temperature 

(0.01  ^  n  ^  1.00;  500°R  *s  T  ^  4000°R)  ,  the  heat  transfer 

w 

coefficient  varies  only  slightly.  As  a  result,  h  may  be 

estimated  without  precise  knowledge  of  r)  and  T^.  However,  the 

wall  heat  transfer  rate  is  a  stronger  function  of  these  two 

parameters.  Curves  of  Nusselt  number  for  heat  transfer  were 

also  presented.  It  is  of  interest  to  plot  Nu  as  a  function  of 

the  local  Reynolds  number  in  the  manner  suggested  by  Reference  17. 

Such  a  plot  is  shown  in  Figure  21  for  several  computed  values  of 

17 

Nu  and  Re.  Knudsen  ana  Katz  point  out  that  the  expression 
which  may  be  derived  for  a  flat  plate  with  constant  wall  tempera¬ 


ture. 


Nu  =  0.0366  (Re)^^^  (Pr)^^^ 


(39) 


has  also  been  found  to  agree  with  experimental  results  for  heat 
transfer  to  a  plate  having  variable  surface  temperature.  Equation 
(39)  is  also  shown  in  Figure  21.  As  may  be  seen,  a  satisfactory 
correlation  is  obtained  when  the  constant  in  Equation  (39)  is 
changed  to  0.0112.  As  a  result,  this  empirical  correlation  may 
be  used  to  obtain  estimates  for  Nu  for  trajectories  other  than 
the  one  considered  here. 
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